Dynamical localization of a particle coupled to a two-level systems thermal reservoir 
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Using the functional-integral method, we investigate the effect of a two-level systems thermal 
reservoir on the single particle dynamics. We find that at low temperatures, within the sub-ohmic 
regime, the particle becomes "dynamically" localized at long times due to an effective potential 
generated by the particle-reservoir interaction. This behavior is different from the one obtained 
for the usual bath of harmonic oscillators and is fundamentally related with the non-Markovian 
character of the dissipative process. 
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I. INTRODUCTION 

The influence of dissipation on quantum tunnehngi 
and on quantum coherence^, has attracted much attention 
during the last decades. In a pioneer work, Caldeira and 
Leggetlji suggested that a bosonic heat bath consisting of 
an inflnite number of harmonic oscillators constitutes an 
universal realization, which can mimic a large variety of 
real environments.— However, a real environment cannot 
always be represented in this way. Indeed, if the parti- 
cle of interest couples to a non linear system, the latter 
may, under certain circumstances, behave as a fermionic 
heat bath. That is the case of a distribution of quartic 
plus quadratic potentials, which within a very well known 
limit ^ effectively acts as a collection of two- level systems 
(TLSs). In such a case the environment truly behaves as 
a spin 1/2 medium. 

A spin-bath composed of an inflnite number of TLSs 
has been mostly considered when the system of interest is 
itself a TLSp providing for instance, a realistic descrip- 
tion of a nanomagnet coupled to a set of surrounding 
spins® and a useful model to describe the loss of quan- 
tum coherence J'- The representation of dissipative envi- 
ronments by a spin-bath has been extended also to situ- 
ations in which the environment is not composed of real 
spins, as e.g. in the description of disordered insulating 
solids- or for explaining the damping of acoustic phonons 
in a nanomechanical resonatorjiS Moreover, it has been 
found that the dissipative dynamics of a single particle 
linearly coupled to a TLSs reservoir is non-Markovian^- 
and that the transport properties strongly differ from the 
usual oscillator thermal bathfi2. Indeed, the optical con- 
ductivity of a set of non interacting particles linearly cou- 
pled to a TLSs reservoir exhibits a remarkable non-Drude 
behavior. In particular, in the sub-ohmic regime the sys- 
tem exhibits a maximum in the incoherent electrical con- 
ductivity at flnite frequency." This kind of behavior has 
been experimentally observed in La2-a;Sra;CuOzpl2, and 
Lai_2;SrMn03jii in situations in which inelastic scatter- 



ing dominates the transport properties. The broad fl- 
nite energy peak observed experimentally suggests a "dy- 
namical" localization of the charged particles, similar to 
the non-Fermi-liquid behavior found in the infrared con- 
ductivity of SrRuOs.'l^' In all the cases discussed above, 
the locaHzation of the particle was attributed to inelastic 
scattering processes because they are strongly enhanced 
as the temperature is raised. This high-temperature 
localization is different from the Anderson localization, 
which tends to be destroyed by inelastic processes. 

In this paper we present a simple model, which leads 
to "dynamical" localization due to inelastic scattering, 
but at low temperatures. At this point, we want to call 
the attention of the reader to our use of the term "dy- 
namical" localization. This term is conventionally used 
in the literature^^ to designate quantum phenomena tak- 
ing place in time-periodic systems, whose corresponding 
classical dynamics displays chaotic diffusion. The phe- 
nomenon described here bears no similarities with the 
latter and for this reason we use the term under quota- 
tion marks. The "dynamical" localization of the charge 
carriers studied here is generated by their coupling to 
a TLSs thermal bath and is in close relation to a non- 
Markovian process at the classical level. At flrst sight, 
a non-Markovian particle dynamics, which implies that 
the particle never reaches equilibrium, seems to be in 
contradiction with localization, which is characteristic of 
insulators. This paper intends also to shed some light 
on this point. In order to achieve our goal, we investi- 
gate the real time effective dynamics of a single particle 
coupled to a TLSs thermal reservoir in the sub-ohmic 
regime. The effective action, which describes the parti- 
cle interacting with the TLSs bath, is obtained using the 
well known Fey nman- Vernon formalism.— 

This paper is divided as follows: In Sec. II we present 
the model describing the particle of interest interacting 
with the TLSs reservoir, as well as a brief sketch of the 
derivation of the effective particle dynamics. In Sec. 
Ill the "dynamical" localization effect is explicitly de- 
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rived and discussed in specific cases within the sub-ohmic 
regime. Finally in Sec. IV we present our conclusions. 



II. THE MODEL 

To begin with, we will describe the particle of inter- 
est coupled to a generic TLSs thermal reservoir by the 
Hamiltonian 



H = 



P 

2M 



u[x) 



N 



/ ^ 2 ^zk - X 2_^Jk<Jxk, (1) 



k=l 



k=l 



where the first two terms stand for a particle under the 
infiuence of an arbitrary potential u, the third term ac- 
counts for the TLSs reservoir, and the last one describes 
the interaction between the particle and the thermal 
bath. Jk denotes the coupling parameter and cr^k/xk 
stand for Pauli matrices. 

The first step is to calculate the reduced density op- 
erator of the particle of interest, which may be obtained 
after tracing out the reservoir degrees of freedom. 



(2) 



The density operator of the total system at time t = 
will be assumed to be decoupled, p(0) — ps{0)pr{Q). The 
reduced density operator ^ can then be written as 

p{x,y,t)= / dx' / dy' psix ,y' ,0)J{x,y,t;x'y' ,0), 



where the super-propagator has the form 

J= r Vx{t') r Vy{t')ei^^°^'=^-'^°^y^'>T[x,y]. (3) 

J x' J y' 

In the expression above, So[x] corresponds to the action 
of a free particle placed in the potential u, while de- 
notes the infiuence functional which describes the infiu- 
ence of the reservoir on the particle dynamics. After 
integrating out the reservoir degrees of freedom and in- 
troducing a set of coordinates corresponding to the par- 
ticle center of mass q — {x + y)/2 and relative coordinate 
^ = x — y,we obtain the super-propagator for the particle 
of interest (see pL3] for details) , 



(4) 



The effective action is given by 



dt' 



Mq{t')i{t')-u(q,i) 



dt"A{t' - t")q{t")£,{t') 



(5) 



where 



A= / dujJ{uj,T)sm[uj{t' -t'% 
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and the functional has the form 

m - /' dt' f dt"^{t' - t"ww), (6) 

Jo Jo 

with 



$ = / dujj{uj, T) cos[cj(i' - t")] coth {huj/2kBT) . 
Jo 

It should be noticed that the kernels A and $ are both 
defined in terms of the spectral density of the thermal 
reservoir™ J{uj,T). 

We now perform an integration by parts to render ex- 
plicit the dependence of the last term in Eq. ^ on the 
velocity. We then find 



5efF 



where 



(7) 



dt' 



t{t' - 1" 



duj^^t^cos[uj(t' -t")]. (8) 

^ 



Although in this form the effective action shows an ex- 
plicit velocity dependent term, characteristic of viscous 
forces, we have obtained also two additional spurious con- 
tributions. The first term in the second line of Eq. ([7]) is 
nothing but a harmonic potential, which can be canceled 
by an appropriate choice of the external potential, 

zi(9,e)-f(0)q^ = f(0)(a;2-2/2)/2. 

This assumption allows us to focus on the dissipative 
effects of the environment. An alternative procedure 
would be to start from a momentum dependent particle- 
reservoir interaction. In this way the effective action 
would immediately exhibit a velocity dependent term, 
without any additional spurious contribution. To con- 
clude the analysis of Eq. |(7| we must point out that 
its last term fluctuates very rapidly for times t :s> l/uj- 
Hence, in principle this term could be neglected within 
the long time regime. However, as we will show below, 
there is no need to introduce approximations because it 
will cancel out naturally when we introduce the initial 
conditions of the problem. The effective action for a sin- 
gle particle coupled to the TLSs reservoir then reads 

SMq, e] = dt' \Mq{t')i{t') - g(0)f (i')C(i') 



dt"r{t' - t")qit")^{t') 



(9) 



Before explicitly solving the equation of motion corre- 
sponding to the action ([9]) , it is convenient to specify the 
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spectral density of the thermal bath in terms of macro- 
scopic parameters. A reasonable assumptions^ for it is 



tanh 



2kBT 



ein-Lo), (10) 



where is a cutoff frequency, 77 is a constant defining 
the coupling strength of the particle to the TLSs, s is 
a number (real and positive) which determines the long 
time properties of the thermal bath, and ujc is some char- 
acteristic frequency introduced in order to make the unit 
of 77 independent of s. Notice that the temperature de- 
pendence in Eq. lfTO|) is crucial for a fermionic heat bath 
because it ensures that the bath degrees of freedom are 
excited as the temperature increases.— In Ref. this 
point was not acknowledged, leading to the wrong con- 
clusion that the decaying term in the particle equation 
of motion is temperature independent. 

Using Eqs. ([9]) and lfTO|) . the classical equations of mo- 
tion for q and ^ can be written as 



27 



27 



dt'T{t - t') q{t') + ^q{0)T{t) 



eE{t) 
M ' 



(11) 



dt'r(t ~ t') e(t') -I- —£,{0)T{-t) = 0, (12) 



where the damping constant is defined as 7 = r]/2M and 
the kernel T is given by 

r(<) = / tanh — — cos{ut). (13) 







2kBT 



Therefore, after tracing out the TLSs reservoir we ob- 
tained an equation of motion for the particle center of 
mass ifTTj) in which the thermal bath has the same effect 
as that of a viscous fluid. 

For any value of s the solution of Eqs. (fTTI) and ifT2|) 
can be written in terms of the kernel Laplace transform 
as 



zq{0) + m 

z2 + 27zr(z)/7r' 



— 2jzT{z)/7r' 



(14) 



(15) 



where 

r(z) 



F n l + s 3+s \ 



tan 



/ hz 



\2kT 



AksTz 



E 



(16) 



In the expression above 2-F1 denotes the hyper- geometric 
function and A„ = (2n — l)TrkBT/h, with n € N. It 
should be notice that the fiuctuating force, given by the 



last term on the LHS in Eq. ifTTj) . does not appear in Eq. 
((HI). This term was exactly canceled by the initial con- 
dition included in the Laplace transform of the damping 
term and therefore there is no need to drop it out by 
assuming the long time approximation. In order to illus- 
trate how the dissipative properties of the TLSs thermal 
reservoir affects the single particle dynamics in the sub- 
ohmic regime (s < 1), lets investigate the simple case in 
which q{0) — and q{0) — vq. 



III. THE "DYNAMICAL" LOCALIZATION 

We start by discussing the s — Q case, in which we 
can proceed analytically a bit further. In this case the 
hyper-geometric function reads 2-^1(1, 1/2, 3/2, —z^) = 
arctanx and the Laplace transform of the damping 
function given by Eq. (flGl) acquires the form 



r = 



AksTz 



tan- (§) 



tan 



-(ft) 



(17) 



In the particular case of 00, the high tempera- 

ture limit of Eq. lfT7|) is r(z) = h/2kBT and therefore 
the effective dynamics of the particle of interest simply 
becomes 



vo 



h-f/2kBT 



(1 



(18) 



This result correctly reproduces the oscillator-bath model 
with an ohmic temperature dependent damping con- 
stant. Indeed, if we assume huj ^ fc^T and fl 00 
in Eq. (fT3l) the damping function becomes T{t — t') ~ 
{hi: /2kBT)5{t — t'), which is an instantaneous func- 
tion, indicating that the dissipative process is completely 
memoryless and the condition of zero frictional force is 
achived only when the particle velocity is zero. Physi- 
cally, this limit corresponds to a weak particle-reservoir 
interaction because most of the TLSs are occupied (on 
average), causing no damping on the particle. Therefore, 
we recover the known result demonstrated in Ref. [l7| . 
namely that when the coupHng between the particle of 
interest and a nonhnear bath is weak enough, the latter 
behaves as a collection of harmonic oscillators. 

However, the features discussed above are not valid for 
all values of temperature and frequency cutoff. In fact, 
we realize that at T = 0, even assuming ^ 00, it is im- 
possible to obtain a damping function without memory. 
In this case the problem becomes non-Markovian and al- 
though the particle becomes localized after some time, 
neither its position nor its velocity ever reach the equi- 
librium. This behavior is illustrated in Fig[T]for finite ft 
and different values of temperature. It is worth to notice 
that in this situation the particle dynamics completely 
differs from its behavior in the oscillator bath model. In 
the later, at zero velocity, the frictional force over the 
particle is also zero and classically, the particle remains 
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Figure 1: Time dependence of the particle center of mass for 
s = and different values of temperature. The continuous line 
corresponds to T = 0.0001, while the dash-dotted, dotted, 
dashed and double-dotted lines are for T — 0.01, T — 0.1, 
T = 0.2 and T = 0.5, respectively. In all cases Q, = 1, ^ — 0.3 
and the initial velocity was taken equal to 1. 



Figure 2: Time dependence of the particle center of mass for 
s = 0.5 and different values of temperature. The continuous 
line corresponds to T = 0.0001, while the dash-dotted, dotted, 
dashed and double-dotted lines are for T = 0.01, T — 0.1, 
T = 0.2 and T = 0.5, respectively. In all cases = 1, 7 = 0.3 
and the initial velocity was taken equal to 1. 



in that state forever. In our non-Markovian situation, 
the frictional force acting on the particle depends on the 
previous velocities with different weights - given by the 
kernel lfT3|) - and the situation of zero frictional force over 
the particle, at a given instant, does not correspond to 
zero velocity. 

It is not difficult to see that at low temperature, the 
kernel oscillates in time, keeping its sign constant and 
therefore the only possibility of getting zero frictional 
force acting on the particle at a given instant occurs when 
the particle changes the momentum direction within the 
time interval. From the physical point of view, this be- 
havior resembles that of a particle confined by a poten- 
tial. Indeed, the particle oscillation around this effective 
potential, which clearly appears for T = 0.0001 in Fig.[Tl 
can be promptly obtained within the long time regime. 
The term cos[uj{t — t')] in Eq. I|13p then oscillates rapidly, 
yielding no contribution to the damping process for long 
times, except when uj {t — t')~^ . In such a situation, 
T{t) becomes finite and time independent, turning the 
damping term into a harmonic localizing potential. As 
the temperature increases, less reservoir states are able 
to play their dissipative role and the particle takes longer 
to get localized far from the origin. 

We can therefore conclude that the particle dynam- 
ics in this situation (s — 0) is completely different from 
the one obtained when the thermal bath is represented 
by the usual oscillator model. Here, memory effects in 
the damping process lead to a "dynamical" localization 
of the particle at a certain distance from the initial po- 
sition, which is proportional to the temperature. The 
strength of the localization potential is determined by 
the ratio of two quantities, namely the thermal and the 
cutoff energies. This quantity measures the number of 
reservoir states which effectively couple to the particle 
and obviously also depends on the total reservoir states 
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Figure 3: Time dependence of the particle center of mass 
for T — 0.001 and different values of s. The continuous line 
corresponds to s = 0, while the dotted and dashed lines are 
for s — 0.5 and s — 1, respectively. In all cases f2 = 1, 7 = 0.3 
and the initial velocity was taken equal to 1. 



determined by the spectral function. 

In order to get more insight on the particle dynamics 
in the sub-ohmic regime we have plotted in Fig. [2]the po- 
sition as a function of time in the specific case of s = 0.5. 
At very low temperatures, the main difference from the 
s — Q case is that the localizing potential strength in the 
long time regime becomes weaker and the dynamical lo- 
calization effect is difficult to observe, see the T = 0.0001 
and T = 0.01 cases, for instance. This is a consequence 
of having decreased the weight of the low energy reser- 
voir states in the spectral function. At the same time, 
as the temperature increases, the number of low energy 
states effectively contributing to the dissipative process 
becomes reduced because several low-energy states are 
occupied (on average) and the particle moves nearly free, 
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see Fig. [2 for T = 0.2 and T = 0.5. This last behav- 
ior differs from the s = situation in which, even at 
high temperatures, there are enough low energy states to 
localize the particle. In general, the particle dynamics 
in the sub-ohmic regime will be described by a function 
which exhibits a behavior in between that of the s = 
and s = 1 Hmiting cases. This behavior is illustrated 
for low temperatures in Fig. [3l The central point in the 
ohmic case (s = 1) is that even at low temperatures, 
there are not enough low energy states that render the 
particle confinement appreciable in the long time regime. 

IV. CONCLUSIONS 

We have studied the real time dynamics of a parti- 
cle coupled to a TLSs thermal reservoir and found that 
within the sub-ohmic regime the particle becomes local- 



ized in the long time limit, oscillating in the real space 
as a consequence of an effective potential generated by 
its interaction with the thermal bath. The oscillatory 
behavior renders the localization "dynamical" and there- 
fore neither the particle position nor its velocity ever 
reach the equihbrium. This behavior is associated with 
the non-Mar kovian character of the dissipative process, 
which in our simple model is provided by inelastic scatter- 
ing of the particle of interest by the TLSs. We hope that 
our findings can be of some help in the understanding of 
transport properties of systems in which the dissipative 
medium seems to be sub-ohmic We also specu- 
late about the extension of the model discussed here to 
a situation in which the fermionic bath excitation has a 
finite gap A(T). In such a case, we expect that the ef- 
fect of "dynamical" localization will start at temperatures 
above A{T)/kB, which is probably more appropriate to 
describe reahstic situations. 
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